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Serai-discrete generalizations of the second order extension of Godunov's 
scheme, known as the MUSCL scheme, are constructed, starting with any three 
point "E" scheme. They are used to approximate scalar conservation laws in 
one space dimension. For convex conservation laws, each member of a wide 
class is proven to be a convergent approximation to the correct physical 
solution. Comparison with another class of high resolution convergent schemes 
is made. 
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Recently, there has been an 

enormous amount of activity related to the construction and analysis of 
"high resolution" schemes approximating hyperbolic systems of conserva- 
tion laws. Some examples of the successful consequences of this 
activity can be found in the proceedings of the latest (sixth) AIAA Com- 
putational Fluid Dynamics Conference, [1], [10], [18]. Extensive 
bibliographies can also be found in these papers. 

Our aim here is merely to construct and prove the convergence of a 
subclass of these schemes approximating scalar convex conservation laws. 
This subclass is based on an idea of van Leer [17]. He christened his 
algorithms MUSCL (Monotonic Upstream Centered Schemes for Conservation 
Laws) schemes, and with mixed emotions, we shall use his acronym here. 


In future work with S. Chakravarthy , we shall extend this construc- 
tion to systems in multi-dimensions, using triangle-based algorithms. 
This work in progress will stress the computational aspects of the algo- 
rithms, especially as they relate to the Euler equations of compressible 
gas dynamics. In earlier work with the same author, [12], we con- 
structed, and proved convergence of, a class of high resolution schemes 
approximating scalar convex conservation laws. We also showed for cer- 
tain high resolution approximations to systems, that limit solutions 
satisfy an entropy inequality. 

We shall consider numerical aproximations to the initial value 
problem for a single conservation law in one space dimension 

( 1 . 1 ) (a) w + f (w) =0, t > 0, -1 < x < 1, 

L» X — 

with a periodic boundary condition: 

(b) w(x + 1,t) - w(x,t) 
and initial condition: 

(c) w(x, 0) = w Q (x). 

It is well-known that solutions of (1.1) may develop discontinui- 
ties in finite time, even when the initial data are smooth. Because of 
this, we seek a weak solution of (1.1), i.e. a bounded measurable func- 
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tion w, such that, for all <6 4 C Q (R x R + ), 

(1.2) (a) SS (wd. + f(w)tf v ) dx dt = 0 

RxR + C X 

(b) lim ! |w(x,t) -w (x)!i . = 0. 

t->0 L 1 

Solutions of (1.2) are not necessarily unique. For physical rea- 
sons, the limit solution of the viscous equation, as viscosity tends to 
zero, is sought. In the scalar case, this solution must satisfy, for 
d 4 C°°(R x R + ), d >_ 0, and all real constants c: 

(1.3) (a) SS ( |w - cid. + sgn(w - c)(f(w) - f(c))d ) dx dt < 0 

u X 

This is equivalent to the statement: 

(1.3) (b) ^ !w -c! + ^ ((f(w) - f(c)) sgn(w - c)) _< 0, 

in the sense of distributions. 

Such solutions are called entropy solutions. Kruzkov has shown in 
[9], that two entropy solutions satisfy 


(1.4) 


iw(x,t 1 ) -u(x,t 1 )!| 1 < | |w(x,t Q ) -u(x,t Q )ii 1 


for all t^ tp. Hence, (1.3) guarantees the uniqueness of solutions 
to the scalar version of (1.2). Existence was also obtained there. 
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Equations (1.3) may be viewed as the statement 


(1.5) ^ V(w) + F(w) £ 0 

for each of the convex "entropy” functions of solutions w to (1.1): 

V(w) = |w - c! 

and their associated flux functions 

F(w) = (f(w) - f(c)) sgn(w - c). 

For scalar convex conservation laws whose solutions lie in the space 
BV, a single entropy inequality (1.5) for any strictly convex V(w) 
and its associated entropy flux F(w), satisfying 

(1.6) F ’ (w) = V'(w)f'(w) 

is enough to imply existence and uniqueness of solutions to (1.1). This 
follows from the results of DiPerna [3]. 

Next we consider a semi-discrete, method of lines, approximation to 
(1.1). We break the interval (-1,1) into subintervals: 

Ij = (x!(j - V 2 ) & £ x < (j + V 2 ) & } 
j = 0, ±1 ±N, with (2N + 1) A = 2. 
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Let x, = j A , be the center of each interval I., with end 
J J 

points *j_ V x J+v 

Define the step function for each t > 0, as 


U^(x,t) = Uj(t), 


for x 4 I.. 

J 


The initial data is discretized via the averaging operator T^, 


T a w q (x) = ^ Xj w Q (s) ds = u j(0) for x 4 Ij, 

J 


For any step function, we define the difference operators 


a ± u j = ±(u j±i - u j > 


D 

+ 


= ± U 


3 


A method of lines, conservation form, discretization of (1.1), is a sys- 
tem of differential equations 


( 1 . 8 ) 




+ D 


■ h j-V 2 


= 0 , 


j = 0, ±1, ..., ±N 


U^(x,0) = T^Wq(x), for x 4 Ij. 
Here, the numerical flux defined by: 
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(1.9) 


h j-V 2 = h(u j+k-l’ * * ' ,u j-k ) * 

for k 1, is a Lipschitz continuous function of its arguments, satis- 
fies the consistency condition: 

h(w,w, . . ,w) = f (w) . 

It is well known that bounded a.e. limits, as A 0, of approxi- 
mate solutions converge to weak solutions of (1.1), i.e. (1.2)(a) is 
satisfied. However, this does not also imply that the limit solutions 
will satisfy any of the entropy conditions (1.5), let alone the general 
condition (1.3). Some restrictions on h are required. 

A simple class of flux functions h, for which (1.8) converges, 
for all f, to the unique entropy solution in L oo (L^(R);[0,T]), as 
/\ 0, for any T > 0, is the class of "E" schemes introduced in 

[11]. Such schemes satisfy the following: 

A consistent scheme whose numerical flux satisfies 


( 1 . 10 ) 


sgn(Uj - Uj^Hh. 


- f(u)] < 


for all u between u and u. is said to be an E scheme. This is 
at present, the most general class of schemes known to converge in the 
nonconvex case. 


It is clear that this class includes the widely known class of 
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three point monotone schemes, i.e. those for which 


h j-v 2 = h( Wi ) ' 

with nonincreasing in its first argument, nondecreasing in its 

second . 


We denote partial derivatives of a numerical flux via 


h(u j+k , ** ,,u j-k+ 1 ) = h r 

Thus a three point scheme with a differentiable flux function is mono- 
tone iff 


Any numerical flux can be written 

(1.10 h j +1 / 2 ■ ValfCuj) * - V 2 Q jt1/2 < 

where Qj + 1/ 2 can be viewed as the viscosity of the scheme [15]. 

One particular three point scheme is due to Godunov, [5], and has 
special significance in this theory. The flux for Godunov's scalar 
scheme can be defined by 

G 

(1.12) h i = h_(u ,u. 1 ) = min f(u), if u. 1 < u. 

2 J J u j_-|£ u < u j J J 



= max f (u) , if u . . > u . 

u _i> u >u . J J * 

J J 

One can thus, [11], characterize E schemes as precisely those for 
which: 


<1 ' 13)(a) h J-V2 - h J-V 2 ’ if u j < u j-1 

(b) lf u j >u j-v 

or, as Tadmor [15] pointed out, those which have at least as much 
viscosity as Godunov's scheme 


(1.14) 


9 j-v 2 i 9 


It follows from [11], Lemma (2,1) 
at most, first order accurate. 


j-V 

, that these approximations are 


f 


Together with an entropy inequality, a key estimate involved in 
many convergence proofs, is a bound on the variation. For any fixed 
t >_ 0, the x variation of U^(x,t) is 


B(y 


= Z 
j 


!A + Uj(t) 


If we can write 


( 1 . 15) (a) 
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then it is easy to show, [12], using an argument of [ 1 9 ] » that, for 




Harten in [7]» pointed out for explicit methods, that this decompo- 
sition could be obtained for schemes which are higher order accurate. In 
our present method of lines context, it involves a five point, con- 
sistent, approximation: 


bu 


(1,17)(a) -St = - A = V/2 - °J-V 2 A.Bj 


with 


b C j+V 2 = C(u j+2 ,U j+1 ,U j ,U j-1 ) - ° 

<0> V/ 2 * D(U J+1' U j’ U j-1’ U j-2 > ^ 0 • 

both Lipschitz continuous functions of their arguments. (See also van 
Leer [16].) 


In addition to (1.16) we have a maximum principle for (1.17), 


(1.18) 


min u. (0) < u .(t) < max u. (0) 
k K J k k 

for each j and all t ^ 0. 

Morover, in [12] we also showed a limit on the possible accuracy of 
approximations of the type (1.17). Any scheme of the type (1.17) is at 
most first order accurate at nonsonic critical points of u. (A sonic 
point u is one such that f ' (u) = 0). Thus, although schemes of this 
type can be made to be as high as third order accurate, Lipschitz con- 
tinuity implies a local degeneracy to first order accuracy at smooth 
maxima and minima. This local degeneracy, together with 3 ome results on 
initial boundary value problems in [6], indicate strongly that overall 
second order accuracy is the best possible. 

Following Harten [7], we call algorithms of the type (1.17), total 
variation diminishing, or TVD schemes. 

It is known, that although (1.17) generates a compact family of 
solutions, U^(x,t), in L 00 (lJ (R) , [0 ,T ] ) , (if w Q -f L 1 (] L 00 f] BV) ), 
limit solutions may not satisfy any of the entropy conditions (1.15), 
and hence may not even be unique (1.12) See e.g. [4], [11], and [15]. 

One way of avoiding this difficulty is to use an inequality obtained in 

[ 11 ]. 


Let V(w) be any convex function. We showed in [11], Section III, 
that, for any solution of any scheme (1.8): 
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(1.19) V( Uj ) + D + F A ( Uj )) = / u J+1 dw V"(w)[h j+1/2 - f(w)], 

J 

where the approximate entropy flux is defined through 

(1 - 20) F A ( Uj ) = F(Uj) + V' ( U j)[hj_i - f( Uj )]. 

Thus, a sufficient condition that any limit solution satisfy (1.5) 
for a fixed convex V, is that 

u . . 

(1 - 21 > X U J+ dw V"(w)[h y - f(w)] < 0 

J 

In order that the above inequality be valid for all convex V, it 
is necessary and sufficient that hj + i correspond to an E scheme, 
which implies again that the approximation be at most first order accu- 
rate. Thus, for the schemes to be constructed in the following section 
we shall only obtain our entropy inequality (1.21) for a single V(w), 
say V(w) = V2W . The following Theorem summarizes the technical 
hypotheses needed for convergence. 

THEOREM 1.1. The sequence of approximate solutions converges a.e. 
— ^ to th e unique solution of the scalar , convex conservation law 

^— *— ) provided that the Initial data is in BV and that inequalities 
(JL‘AZ.) a °d (1.21) are valid for £ single convex V(w). 

The proof is analogous to that of Theorem (4.1) in [12], 
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IX. Construction of TVD MUSCL Schemes. The MUSCL scheme as dis- 
cussed e.g., in [2] and [8] is a second order accurate extension of 
Godunov's method [5] that is based on ideas expressed by van Leer in 
[16]. 


We assume that u.(t) is known. The first step in MUSCL is the 
<3 

reconstruction of a piecewise linear description of the solution. In 
the interval 1^ the result of this operation is 


(2.1) z(x,t) = Uj(t) + (x - Xj)Sj(t), for (x - x^) < A /2. 


Here s.(t) 
J 


are slopes that satisfy 


(2.2) 3^ = w x (Xj,t) + 0( A ) 

subject to some monotonicity constraints discussed below. 

Godunov's scheme, in our semi -discrete context, is the following. 
The numerical flux h G (Uj +1 ,u^) is computed by solving the Riemann 
problem, i.e. the initial value problem (1.1) with initial data 


w s Uj(t), for x < x j+1/2 
w-u j+1 (t), for x > x j+1/2 - 

The resulting unique entropy condition satisfying solution is a 
function of the type 
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A closed form for w was recently obtained in [11]. 


Then Godunov's flux is defined through 

(2.3) h^ = f (w(0) ) 

which is the same as (1.12). 

A fully discrete, explicit in time, Godunov scheme has the same 
numerical flux, after we impose a CFL restriction which prevents the 
interaction of solutions of adjacent Riemann problems The literal 
second order accurate extension to this fully discrete approach would be 
to compute the solution to the initial value problem (1.1), with initial 
data 


w r u.(t) + (x - x.)s.(t) , for x < x . i. 

J J J j+'/2 

8 U j + 1 (t) + (x " x j )s j + 1 (t) ’ for x > x j+ 1/ 2 - 

Find f(w(x j + i^»s)) for t + X /\> s >_ t, with X > 0 satisfying a CFL 
restriction, then compute 

(2 - M) T7\ *t +> ^ f<w(x j + i/ 2 ' 3 » ds - 

This is the MUSCL numerical flux. 



Unfortunately, obtaining the exact solution to this nonlinear ini- 
tial value problem with piecewise linear initial data is a non-trivial 
business, even in the scalar case. However, at s = t + , x = x j + 1 / 2 » the 
solution is the same as for the Riemann problem with initial data: 


w r u j(t) +Y s j( t )» for x < Xj + i/ 2 

“ ?u j*i for 

which is easily calculated, e.g., from the formula in [11]. Thus, the 
semi -discrete MUSCL extension of Godunov's algorithm comes through 


c5u 


(2.5) 


'St = ~k ^- hG(u j + i ■ ^ S j+r u j 


We generalize this still further as follows. Let 

h(u. . ,u.) = h. i. be an arbitrary first order accurate numerical flux 
j+1 j J+'/2 

function. Then our generalized MUSCL algorithm is merely: 


( 2 . 6 ) 


c)u . 

■St 1 


-k ^- h( V i *z V- 


We have: 


LEMMA 2. 1. At points w(Xj,t) = Uj(t) in £ neighborhood of which 

h(u. .,u.) is C 2 with Lipschitz continuous partial derivatives , and 
* J 

A3 i+ 1 /13. 

-rr T 7 ~~ = 1 + 0(A) = 7 ? — the algorithm (2.6) is at least second order 
j Uj -- 
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accurate for smooth functions w 


Proof. Let 


(2.7) 


u j.v 2 = V 2<Vl * U j>' 


We shall show: 


(2.8) A_ (h(u j+1 -^s j+1 , Uj Sj) - f(u j+1/2 )) = o(A 2 ) 


at these points. 


Using consistency of h(Uj +1 ,Uj), the left side of (2.8) can be 
written as 


/\ 


Jo h , 


A* u j 

u j-.v 2 ♦ y -r^ 


As. ’ 


^ u j 


lU j *2 S i 


dy 




A s 




+ h o 


u < + i, + y 


+ 2 j’ j+V 2 

The result is immediate. 


A, u [As. )) fA u f A s i *j 

2 JJ l 2 M 


Next we make an observation about the viscosity of this scheme, 
Qj+1/ 2 * as defined in (1.11). 

LEMMA 2.2. I£ h j^s monotone , then the viscosity of algorithm (2.6) 


is a decreasing function of and J 

57^— 57*7 



The proof comes directly from (1.11) 


Thus 


s . s . . 

, as expected, if -r — “ — and 7 — ** — 

b. Uj a u j 


are both restricted to be 


positive, the most viscous, hence the least accurate, case occurs when 
these values are both zero. The scheme then degenerates to first order 
accuracy. We can compress or smear out the solution locally by increas- 
ing or decreasing these ratios. 

We wish this scheme to be TVD, which will give a natural restric- 
tions on these above ratios. 


LEMMA 2.3. If h i£ £ flux corresponding to an E scheme , then the 
scheme is TVD if 


(2.9) 

for each j . 


0 < 


As, 


~ u ; 


^ 3 j+1 


< 1 


Proof. We have 


(2.10) -f Sj_,) 


s a [hl Vi 


A* 


S J> - f(u j *^ S J> 


As, 


. f(Uj .£ Sj > -h(Uj -f Sj)] 
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A Ch(u j - V’ u j +/ i 3 


) - f (u - 


/\ s. 


* f(u j --t 1 ’ - h(u j -^-vi 

Thus, in (1.15) we can write: 


(2. 1 1 ) (a) C 


g 

J +1 /2 = " E Uj) [(u j+1 " “2~ ^ u j + ^ s j ) “ f(u j + § 3 j ) 


f( u. + ^ 3j ) -h(u. -£s jfUj + ^ Sj 


)] 


(b) 


^ g 

D 3-V 2 = ^ »j) tMu j ' T^J *£v - f(u j -f 3 j> 


* f(U j -2 s j ) - h(u J - 2 3 j>Vl *2 3 J-1 


)] 


Both C j + i /2 » D j_1/ 2 are nonnegative because of (1.10) and (2.9). 

The restrictions (2.9) can be relaxed somewhat if h is monotone; 
we have: 

LEMMA 2.4. I£ h ijs a flux corresponding to a monotone scheme , then 
the scheme is TVD if 


(2. 12) (a) 1 > (a J+1 ~ a i ) A 


2 7 ^ Uj near P°ints where h 1 j! 0 
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( — S . . + S , ) A 

1 :> 2 A ' u near p0 * nts where h 0 M 

^ j 


Proof. We have 


(2.13) i/VhOij -f s j-Vl 


. As u , As^ As. As 

= A [ »<Vi - -2-’ u j * -r 4 - h(u j - -r^j ♦ -r 1 ’ 1 


As, As, 


As, As. 


+ A Ch(u j - -r 1 *^ + -r 20 - h(u j - + ~ 2^- )]< 


Thus, in (1.15), we can write: 


(2. 14) (a) C >1/2 ^ - A(^ ~ u -) h 1 (g j' u j - 2< 3 j.l ‘ 3 j )] 


As, 


<b) V/ 2 * S’sbv h ° lU j ' “j * 2 <3 J ' 3 j-1 


Both Cj y and y^ are nonnegative because of monotonicity 


and (2.12). 


Remark 2.5. These last two Lemmas can be made "local," i.e. in 
regions of monotonicity we may relax the restrictions from (2.9) to 
( 2 . 12 ). 
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.-j - J-j J S c he me s . In order to prove convergence 

for convex f, we need only verify the discrete entropy inequality 
(1.21) for the schemes constructed in the previous section (Theorem 
(1.1)). We shall do this for the entropy 


V(w) = V 2 w 2 

and for schemes satisfying the TVD hypothesis of Lemma (2.3), in addi- 
tion to some other restrictions near sonic shock points. 

Suppose Uj < Uj +1 , i.e. we are at a rarefaction. We then have, 
letting u j + l /2 = V 2 (u^ + u j+i^ as ab c» ve 


u -ui &s. , As. 

(3.D S u J dw[h(u j+1 - — ! -,Uj + -yJ-) - f(w)] 

J 


3 -i + 1 3 

= <&. - -j—. “j . - r(u j + 1^)) 


* A.J* 1 d “ cr(u j.i /2 > - 


= [I] + [II]. 

Now, using the hypotheses of Lemma (2.3), we have: 


(3.2) 


[I] < 0 
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Next we have 


( 3 . 3 ) [II] = ^ j+1 dw f ” j+1/2 ^j(s - w)J f'(s) 
J 


■fuj * 1 d ” ( V /2 - "> f ' <U j» 1 / 2 > 


S dw / J + ^2 ( 3 _ w ) f”(s) ds 
v u. *w 

J 


- J* j+1 dw / j +1/ 2 (s - w) f"(s) ds 


= - \ J* U J+ ds f"(s) (minis - u j+1 l,|s - Uj j T < 0. 

Suppose u. > u. i.e. we are at a shock. Let Q . i, be chosen 
J j+i J+/? 


so that 


f'lwKw - a jt y 2 ) d« = 0 


< 3 .«> 


f.(.) d» 

V/ 2 = ~fa 57^ 

i u J+1 f'(w) ^ J 

J 
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where F(u) is the entropy flux corresponding to our entropy, 

V(u) = V2U . Of course we need the denominator in (3.4) to be nonvan- 
ishing. In fact we restrict ourselves to the situation where 

<3 - 5> Vi * V* 1 V 

We note 


(3.6) 




J+1 


(w " “j+l/^ f,(w) dw 
f(uj) 


' u f’ <“ - V/ 2 > - f'fvV 1 d “ 

Thus, we have 


(3.7) 


sup f"(u) 

1 5 U. ..<U<U, 

’V / 2 - W 2 ! < J V • 


A- f(U 1 ) 

where s, 1r = — . 

j+V 2 u j 


We further restrict the so that 


(3.8) 0 < -/iSj < 2 max(min((u J - d J+1/fe ) . (Dj.-i^ - Uj)),0). 

In view of (3.7), this restriction does not affect the second order 
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accuracy of the scheme, except at sonic points 


We now write 


(3.9)/ u j+1 dw[h 

j 


A 3 ,1+1 

U j+1 " 2 * j + 


As, 


j 2 


- f(w)] 


= ^ u jt h 


„ ^JLi 

U j+1 " 2 ’ U j + 2 


“ f( W 2 )] 


+ sj dw crew) - f(a j+1 )] 


j+i 


= [III] + [IV]. 

Now, the assumptions (3.5), (3.8) imply 

(3. 10) [III] < 0. 

An argument analagous to that of (3.3) suffices to show that 
[IV] < 0. 

Finally, suppose u^ +1 < u^ and (3.5) is violated. The inequali- 
ties of Lemma (2.1), and (3.8) imply that s^ and Sj +1 vanish. Then 
inequality (1.21) is immediate, since h is an E flux. 

We may summarize all this in 
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THEOREM 3.1. The sequence of approximate solutions satisfying 
(£.6) converges £.e. to the unique solution of the scalar convex conser- 
vation law 0.1) provided that the initial data is in B.V. and that for 
each j 


0 < =■*- 


> /■ \ 

-^ U j" 


< 1 


0 < 


& 3 , 


U 1 


< 1 


and, in addition if u, > u. ... 
j j+r 


- /\ Sj < 2 max (min ((u. - Qj^), - u^)) t °)» 

JV. SteadyDi^reteShocks . We now check for the existence of 
discrete, steady, shock solutions to the convergent MUSCL schemes con- 
structed in the previous section, based on the Engquist-Osher flux, [4], 
Since the scheme satisfies a maximum principle (Lemma (2.1) of [12]), 
any profile must be monotone. 

Let f"(u) > 0, and for simplicity, we take f ' (0) = 0 = f(0). In 
thi3 case, the monotone E-0 scheme becomes 


(4. 1 ) (a ) h E0 (u j+1 ,u j ) = f_(u j+1 ) + f + ( Uj ) 
where 
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(b) f Ju) r f(u) if u < 0 

f_(u) = 0 if u > 0 

(c) f + (u) - f(u) if u > 0 

f + (u) - 0 if u < 0 

L R 

Let u , u be the left and right states for a physically correct, 
steady shock, i.e. 

f(u L ) = f(u R ), u L > 0 > u R . 

A steady discrete shock {u j } will satisfy 

(4.2) (a) lim u . = u L 
J-^-oo J 

(b) lim u . = u R 
j-*M-00 ^ 



for all j. 
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We shall seek steady, discrete, shocks of the same general form 
obtained in [4] for the E-0 scheme, c.f. also [12], section V. 

These are 


(4.3)(a) Uj - u L , j < -1 

(b) u j r u R , j > 2 

(c) u L > u 0 > 0 > u 1 > u R . 

For the first order scheme, can be viewed as a smooth function 

of satisfying the above, and, in addition 


f(u Q ) + f( Ul ) = f(u L ) 


For the present scheme, we shall get a different one parameter family of 
intermediate states. 

It follows for all j * 0, that (4.3) implies (4.2) (c) . 

For j = 0, we have the following equation: 


t 


«. / 


s 


(4.4) 0 = f 


A s 1 


U 1 ~ 


+ f 


s t 


u 0 + 


- f(u u ) = G( Ul ,u 0 ) 


One special solution is, again, u 


1 


For (Uq,u^) close to (0,u K ) 


t 



and satisfying (4.3)(c), we have, 
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from the hypotheses of Theorem (3.D: 


(4.5) s, = a 


Cu - u,) 


i = a v 2 A 


s Q = 0 


Here 0 < ai. is assumed, in addition, to be a C function of 

" 2 * 

R ba Vp 

u 1 near u 1 = u , with = 0. 


A straightforward calculation gives us 


(4. 6 ) (a) (u R ,0) = f.(u B > + -^J < 0 


(b) (u R ,0) = 0 


(c) (u B ,0) 




«> 5 ^ * 0 


X 2 r 

(e) “ = f..( 0 ) > o 

bu o 


A simple application of the implicit function theorem gives us: 
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THEOREM 4.1. Given the hypotheses of Theorem (3.1> and the mild 
technical assumption concerning a^, there ex ists a family of sharp , 
discrete, shock solutions to (2.6) of the t^ge (4. 3) , with Ul a smooth 
f unction of Uq, and 0 < u^, small enough . By symmetry , the same is 
true, with u Q a smooth function of u 1 < 0, and -u 1 small enough . 


~ In 

[12] we constructed a class of convergent high resolution schemes, using 
flux limiters and second order accurate upwinding. Our first step was 
to use a first order E scheme flux, h(u j+1 , Uj ), to construct a second 
order accurate TVD scheme. The resulting scheme is 


(5.1) 

with 



A 


/X.H(u 


j+2, U j+1’ u j 


’ U j-1 


), 


(5.2) H(u j+2 , Uj+1 , Uj , Uj-1 ) = h(u j+1 , Uj ) - V 2 * (Rj + 1 )(h(u j+1 , U j) - f( Uj ) 


+ V 2 i>(Rj)(f(Uj +1 ) -h( U j +1 , U j)) 


Here 


(5.3) 


,♦ f(u .i> - h(u i- u i-i> 
J - h(u j+1> 


tl<u .i.i* u i ) - »“i> 

3 h(Uj, u J _ 1 ) - ftuj.,) 
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The slope limiter i) (R) satisfies i>( 1 > = 1 and the following 


inequalities 


(5.4) (a) 


1 + V 2 


CRT) 


> 0 


(b) 


i + v 2 


? (R ? 


R 


- * CR* >! 


> 0 


Various slope limiters have been developed. See Sweby [13]* for a 
numerical and theoretical analysis of their properties. Perhaps the 
simplest example is the so called min-mod, [14]: 


(5.5) 


(R ) = max (0,min(R, 1 ) ) . 


We have the following analogue of Lemma (2.2): 

LEMMA 5.1. If h is an E flux , then the viscosity of (5.J.) is a 
decreasing function of »>. 

Next we compare the schemes (2.6) and (5.1). If f'(w) t 0 and w 
is smooth at x^, then it is easy to see that 


H(u j+ 2 ’ Vi' V Vi’ * h |Vl ‘ 2 * • U J * ~2 


2 2 
+ o((^_ Uj) + u j) ) 


if 
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(5.6) Asj = i)(R") A. u j + °<(A. u j) 2 + u j )2) ' lf ^(w.w) / 0 

and/or 

Asj = 'XR*) /\. u j + 0((^ Uj) 2 + (^ Uj) 2 ), if h Q (w,w) * 0 
Dropping the quadratic terms, inequality (2.9) gives 


(5.7) (a) 


0( RT) 

°< * ^(RT) £ 1 if h,(w,w) / 0 

r: j 1 


(b) 


<>(R*) 

0 < — , ?(R*) < 1 if h n (w,w) 4 0 

rT J u 


These inequalities are a bit more restrictive than those needed in 
the monotone case. There we have 


(5. 8) (a) 1 > 1/ 2 


*(RT) 

j 


if h.j(w,w) £ 0 


(b) 


1 > V 2 


'MR*) 

— f- + »>(Rt ,) 


‘j 




if h Q (w,w) ^ 0 


which are the same as inequalities (5. 4) (a) , (b) . 


For the entropy condition to be proven valid, the flux was modi- 
fied. We replaced H by H ac (where a.c. stands for artificial 
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compression), as follows 


First we modify the flux differences in (5.2): 


(5.9) (a) (h(u^ +1 , Uj) - f(Uj))‘ 


= (h(U Jt1 , Uj) - f (Uj)) 


h,(u r u,))^ U,) 
hlu jt1’ u j> - f(u j ) 


(b) (f(Uj +1 ) - h(Uj +1 , Uj)) 


(f(u. + 1 ) -h(u j+1 , Uj)) 


a^ttyw y>/y u. 
f ( Vi> - h( Vi’ u j> 


Here a 7 i, are both positive, chosen first so that the quantities 
J +l /2 

in the brackets in (5.9) are between 0 and 2, (inequalities (5.13) 
below) . 


We next let: 



(f(u ) - h(u . 
(f(Uj +1 ) - h(u 


. u 


j+1 



, Uj )) 


M 


u,) -f(u t )) M 
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Then we define our scheme: 


( 5 . 10 ) 


Tx~- H “ <u j.a’ Vi- V Vi> 


with 


(5.H) H aC (Uj + 2* « J+1 , U jf Uj.,) 

H M 
= h(u j+1 , Uj) - V 2 (h(u j+1 , Uj ) - f(Uj)r 

M 

+ v 2 W+) Cf<u j+1 ) - h(u j+r Uj )) M . 

The resulting scheme is convergent if the following inequalities 
are valid: 


u . . 

(5.12) X J+ f"(w)dw ((V 2 ( A u . ) ) 2 - (w - V 2 (u. , + u.)) 2 ) 
j * J J+ ' J 


" V'/ 2 'A u j> th i ( Vi- u >i> - W “j* 1 


‘ 3 j+i/2 <4. u j >2 th o ( Vi'Vi) ' h o (u j’ V 1 - 0 


, _ (h(u , u .) - f(u .) ) 

(5. 13) (a) ia. +1/2 (hi(u j+ i, u J+1 ) - h^u., Uj ))| < - — ^ - 1^3 i_ 


(b) 


ia j+V 2 (h 0 (u j+1* u j+i ) _h o (u j’ u j ))! i 


f(u > i ) -h(u i+1> Ul ) 

^u. 
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and if i$(R) is defined by (5.5). 


Examples are given in [12]. 

It is interesting that this case allows compression to be added at 
shocks, i.e. the restrictions on the a"^ allow negative viscosity to 
be added to (5.2) if > Uj +1 , while, again for the proof of conver- 
gence, positive viscosity must be added for rarefactions, Uj < Uj + t In 
the present MUSCL case the situation is (unnaturally) reversed. A com- 
parison of the convergent MUSCL scheme with this one yields two main 
points: 

(1) Away from sonic points, the fluxes differ only in 0( (A t u^) ) 

terms 

p 

(2) At sonic points, the fluxes differ in 0((^ Uj) ) terms. 
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